Joint quantum nondemolition measurements of qubits: beyond the mean-field theory 
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We propose an approach to nondestructively detect A'^ qubits by measuring the transmissions of a dispersively- 
coupled cavity. By taking into account all the cavity-qubits quantum correlations (i.e., beyond the usual coarse- 
grained/mean-field approximations), it is revealed that for an unknown normalized A'^-qubit state JV'iv) = 

X^j,^o Pk\k)N, each detected peak in the cavity transmitted spectra marks one of the basis states |fc)jv and 
the relative height of such a peak is related to the corresponding superposed-probability |/3fc p. Our results are 
able to unambiguously account for the intriguing multi-peak structures of the spectra observed in a very re- 
cent circuit-quantum-electrodynamics experiment [Phys. Rev. A 81, 062325 (2010)] with two superconducting 
qubits. 
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PACS numbers: 42.50.Hz, 03.67.Lx, 85.25.Cp 



Introduction. — It is well-known that the readout of qubits 
is one of necessary steps in quantum information processing. 
Phenomenally, the information stored in an unknown N-qubit 

2"— 1 
quantum state I -0 at) = X)a:=o /3a: | A;) w can be partly extracted 

by directly applying the standard von Neumann projective op- 
eration P = ^k |^)Afjv(^| to the quantum register (111- After 
such a projection, the register is collapsed to one of the com- 
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7=1 \'^j)N,ct = 0, 1} with a probability |/3^p. This is a 
directly projective measurement (DPM) and the register is de- 
tected. Typically, DPM has been utilized to extract the bi- 
nary quantum information stored in |2], such as trapped ions. 
Cooper-pair boxes, and the cuiTent-biased Josephson junc- 
tions, etc.. Essentially, due to the inevitable back actions of 
the the measuring apparatus, the detected 1/3^^ is always less 
than its expectable value |/3fcp. This means that the DPM is 
not an ideal method to extract the quantum information in an 
unknown quantum state i^ ■ 






FIG. 1: (Color online) QND measurements of A'^ qubits (with two 
levels: {|0j), |lj)}, j — 1,...,N) by measuring the steady-state 
transmitted spectra Sff{LJL) of the dispersively-coupled cavity (with 
frequency cjf and decay rate k) driven by a frequency-controllable 
external field E{ujl)- 



Alternatively, indirectly projective measurements (IPMs) 



can also be utilized to achieve the measurement of the qubits, 
where another coupled system instead the qubits-selves is de- 
tected. Typical examples include, e.g., dc-SQUIDs for the 
inductively-connected Josephson flux (persistent) qubits |Q], 
optical cavities for the containing atomic qubits Q, and 
Cooper-pair box for the nanomechanical resonators |6], etc.. 
A remarkable advantage in the IPM is that the back actions 
from the detected system could be minimized. If the con- 
dition [iJjv,i//] = is satisfied (which means that the dis- 
turbance of the detector D on the qubits is negligible), then 
the relevant IPM further becomes a quantum nondemolition 
(QND) measurement fr]. Here, Hn is the Hamiltonian of the 
N-qubit register and Hi the interaction between it and the de- 
tector. Note that the term nondemolition does not imply that 
the wave function of the register fails to collapse due to the 
measurement [7]. In fact, for an unknown input state \iIin), 
after the QND measurement the N-qubit register will be au- 
tomatically collapsed to one of its computational basis \k)p4 
with an ideal probability |/3fep. Thus, QND measurement is 
a conceptually ideal projective measurement; the successive 
QND measurements on the same register will give the same 
result. 

As it can be easily detected with the cuiTent technique, 
driven cavity has been widely utilized to achieve the desirable 
IPM of the dispersively-coupling qubits. Experimentally, it is 
not difficult to probe the resonance frequency of a driven cav- 
ity by detecting the transmitted signals. If the qubits are dis- 
persively coupled to the cavity mode, then the cavity is pulled 
by the qubits, depending on the state of the qubits. As a con- 
sequence, by detecting the shift of the central frequency of 
the driven cavity mode the QND measurement of the qubits 
can be achieved. This idea has been experimentally demon- 
strated by the cavity QED experiments with few qubits, and 
single basis states (i.e., computational basis) of the atomic 
and superconducting qubits had been experimentally distin- 



guished fs, 9] . Next challenge is to completely characterize 
an unknown N-qubit state IV'at) by nondestructively measur- 
ing an arbitrary superposition of the single basis states with 
ideal precisions. 

Considering the practically-existing dissipation of the de- 
tector (i.e., cavity) and also the statistical quantum correla- 
tions between it and the iV— qubit quantum register, in this 
letter we show that an unknown quantum state | V'jv) can be ef- 
fectively nondestructively detected by the realistic QND mea- 
surements. Our proposal still works for the mixed states, 
and thus could be utilized to explain the detected multi-peak 
structure in the most recent circuit-quantum-electrodynamics 
(circuit-QED) experiment III2II . where the detected qubits 
could be prepared/decayed at various superpositions of all the 
possible basis states. 

Generic model. — The system proposed to nondestructively 
detect an iV-qubit state is schematized in Fig. 1, wherein N 
non-interacting qubits are dispersively coupled to a driven 
cavity. Certainly, the preparation of the initial state of the 
qubits and detection of the driven cavity are repeatable. With- 
out loss of generality, we assume that, the coupling strength 
gj and the detuning A^ = \u>f — ojj \ between the jth qubit 
and the cavity, and the detuning Aij — \uji — loj \ between the 
jth and jth qubits satisfy the condition 
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«1, ^^J = 1,2,...,7V. (1) 



This is to realize the dispersive interactions between the qubits 
and cavity, and to assure that the ith and jth qubits are decou- 
pled effectively from each other Also, the decay rates {7^} of 
the qubits should be significantly less than that of the cavity, k, 
to ensure that the detected state of the qubits has sufficiently- 
long lifetime. In fact, all the conditions listed above are prac- 
tically satisfied in the current typical circuit-QED systems |9]. 
Under the rotating-wave approximation and in a framework 
rotating at a frequency lol, the process for nondestructively 
measuring the unknown A^-qubit state could be described by 
the following master equation 



Pn = -i[HN,PN] + -{2apNo) - a^ap^ - pNa^a),(2) 
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Hn = Sa^'a + -;^'^^]<^j - o)a. ^ rjcr| + e{a) + a), 
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with Vj ~ 9]! ^j and e being the effective strength of the 
driving. Also, a)j = loj — Y j is the renormalized transition 
frequency of the jth qubit, and (5 = w/ — cj^ is the detun- 
ing between the driving field and the cavity. For the ideal 
readouts, the dissipation of the qubits is assumed to be neg- 
ligible. As a consequence, the expectable values of all the 
related qubit-operators, i.e., a^;, ffo'fe (J 7^ fc)> and the iV- 
body ones Hi^i '^f are kept unchanged. Our central task is 
to calculate the steady-state transmitted strength S^{ujl) of 
the driven cavity. This quantity is essentially proportionate 
to the number of the steady-state photons in the cavity, i.e., 
Sff — (ata)^/e^, which is determined by the following dy- 



namical equation 



dt 



~ ~K {a^a)N — 2elm(a) 
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Here, {a)]y = Tr(apjv) is further determined by 
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Neglecting all the statistical quantum correlations between 
the cavity and qubits, i.e., under the usual coarse-grained (or 
mean-field) approximation (CGA), see, e.g., lUOll . we simply 
have (CT|a)jv ~ (CT|(0))7v(a)7v- Then, by finding the steady- 
state solutions to the Eqs. (3-4), one can easily obtain an ap- 
proximate transmitted spectrum: 



S-N'i^L) = |k - K - ACo^)]' + (f )'} 



(5) 
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with Acjjv = Si=i rj(cr|(0))Ar. This indicates that, com- 
pared to the spectrum for the empty cavity (EMC) transmis- 
sion, the qubits only shift the central frequency with a quan- 
tity Acjjv and the single-peak shape is unchanged. However, 
the above CGA is unnecessary and the two-body cavity-qubits 
correlation functions ((T^a)Ar can be further determined by 
solving the following dynamical equation 
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d{cr^a) 

dT 

(6) 
Note that the three-body cavity-qubits correlations {a^afa)M, 
introduced above, is related further to the four-body cavity- 
qubits correlations: {a^ af a^d) n , m = 1,2, ..., N, etc.. Gen- 
erally, the fc-body cavity-qubits correlations are related further 
to the (fc + l)-body cavity-qubits correlations (i.e., k qubits 
correlate simultaneously to the cavity), and thus a series of 
dynamical equations for these correlations will be induced. 
Fortunately, due to the fact that afa^ = 1 for I — m, these 
equation-chains will be automatically cut off and ended at the 
{N + l)-body cavity-qubits correlations. Then, all the inter- 
ested statistical quantum correlations in these equations can 
be exactly calculated, and consequently the transmitted spec- 
tra can be obtained beyond the usual CGAs. It is emphasized 
that the spectral distribution 5^(wl) including all the cavity- 
qubits quantum correlations may reveal 2^ peaks for an iV- 
qubit state superposed by 2^ basis states. If the detected state 
is just one of the basis states (not their superposition), then 
Sf^{ujL) reduces to S^{ll!l) (with a single-peak structure) 
and the cavity-qubits quantum correlation vanishes. 

Demonstrations with experimentally-existing circuit-QED 
systems. — Our generic proposal derived above could be 
specifically demonstrated with various experimental cavity- 
qubits systems, typically the circuit-QED one 113, 1 1411 . In 
this system, the cavity is formed by a coplanar waveguide 
(of the length at the order of millimeters) and the qubits are 
generated by the Cooper-pair boxes (CPBs) with controllable 
Josephson energies. At a sufficiently low temperature (e.g.. 



N- 



< 20 mK), the coplanar waveguide works as an ideal super- 
conducting transmission line resonator (i.e., cavity). Exper- 
imentally llioll . the decay rate (e.g., k ^ 2tt x 1.69MHz) 
of the cavity is about ten times larger than that of the CBP- 
qubit (e.g., 7 = 27r x 0.19MHz) UIl]. Also, by adjust- 
ing the external biases, the CPB-qubits could be either cou- 
pled to or decoupled from the resonator, and the required ini- 
tial state preparation and detection can be robustly repeated. 
For the EMC case, the steady-state solutions to Eqs. (3-4) 
can be easily obtained and the transmission spectrum reads 
SI,'{ijJl) = [(wi - ujff + (K/2)fJ-i. Obviously, this is a 
well-known Lorentzian lineshape ||9[] centered at ujf, with the 
half-width k. 

For one qubit case with A^ = 1, the steady-state transmis- 
sion spectrum of the cavity is expressed as 



sriu^L) 



{ujL - LOff - 2{lol - ^/)riZf ^ (0) + Ai 



with Ai = r? + (k/2)2 and zf ^(0) = Tr{pi(0)crf} == 
2|/3ip — 1, for the unknown qubit state |^i). This is ev- 
idently different from the SI^Ilul) derived under the usual 
CGA. Obviously, the spectrum function ^^(wl) predicates 
that two transmitted peaks could be found in the spectrum. 
Specifically, using the parameters in the experimental circuit- 
QED with one CPB-qubit [lH: {ujf,uJo,K,g) = 2tt x 
(6444.2, 4009, 1.69, 134)MHz, we plot respectively the spec- 
tra Sf^{uJL) and Sf^{ujL) in Figs. 2(a) and 2(b) typically for 
|/3;^|2 = 0,0.25,0.5,0.75, and 1, respectively. For contrasts, 
the spectrum of the empty cavity (black line) is also plotted 
in the figures. We make two remarks. Firstly, if the qubit 
is at one of their basis states (i.e., either |0)i or |l)i), then 
Sf^iujL) and S'f'*(aji) give the same single-peak distribution, 
which has been experimentally demonstrated [9]. Secondly, if 
the qubit is prepared beforehand at the superposition of its two 
basis states, then 5f'*(aji) shows still the single-peak struc- 
ture (when |/3i| = 0.5, S'f*(ajL) superposes the Sq''{ujl)) but 
Sf''{ujL) predicts two peaks: the locations of the central fre- 
quencies are unchanged, but their relative heights equal re- 
spectively to the superposed probabilities of the two basis 
states. Therefore, 5f*(aji) (rather than ^^^(wl)) provides 
the messages of all the diagonal elements of the density ma- 
trix pi. These predictions should be easily verified with the 
current experimental technique, once the qubit is input at an 
arbitrarily-selected superposition state. 

Similarly, for N = 2 case lll2l - ll4ll the steady-state trans- 
mitted spectrum can still be analytically obtained: 
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FIG. 2: (Color online) Left: Spectral distributions Si"{ujl) (a) and 
Si''{ijJl) (b) of the cavity with a single qubit prepared at the state 
\ip)i =/3o|0)i+/3i|l)i,typicallyfor|/?i|2 = 0, 0.25, 0.5, 0.75, and 
1, respectively. Right: Spectral distributions S2°(aj_t) for the two- 
qubit register prepared at IV'2) with (|/3o|M/3i|M/32lM/33p) = 
(1, 0, 0, 0), (0.34, 0.66, 0, 0), (0.47, 0, 0.53, 0), (0.2, 0.2, 0.26, 0.34), 
respectively. In contrast, transmission spectrum of the empty cavity 
(EMC) So''{ujl) is also shown by the black line. 



(7) Tr{p2(0)(Tfa|}, zf\o) = Tr{p2(0)a|} and D 
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1,2, respectively. It is seen that, the spectral dis- 
tribution S'|*(wl) may reveal four peaks, but 51* (a;^) 
always shows one-peak structure. The current circuit-QED 
experiments with two CPB-qubits III2II could be utilized to 
verify the multi-peak spectral distributions predicted above, 
once the CPB-qubits are prepared at the superpositions 
of their basis states. Specifically, with the experimental 
parameters |[Il[ll: (w/ = w^Fi = x^,r2 = X^,k) = 
27r X (6.806, 0.013, 0.004, 0.001)GHZ, one can plot the 
spectral function 5*1* (wl) in Fig. 2(c) for the selected two- 



qubit state \^J2) with: (|/3o| 



M/32|M/?3n = (1,0,0,0) 
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(black line), (0.34,0.66,0,0) (blue line), (0.47,0,0.53,0) 
(red line), and (0.2,0.2,0.26,0.34) (pink line), respectively. 
From these numerical results one can see that, if the two-qubit 
is prepared at one of the basis states, e.g., |00) here, then the 
transmitted spectrum of the driven cavity shows a single peak. 
While, if the two-qubit are prepared at the any superposition 
of their basis states, then the detected spectra should reveal 
multi-peak structures, i.e., two peaks for black and red lines 
and four peaks for the pink line. 

Generally, if the N-qubit register is prepared at the super- 
position of A/(< 2^ — 1) basis states, then the cavity could 
be pulled by M forms and thus there are M possible shifts 
of the cavity resonance frequency. As a consequence, the de- 
tected cavity transmitted spectrum will reveal M peaks; the 
superposed probability of one of the basis state determines 
the weight for pulling the cavity and thus the relative height 
of the coiTesponding transmitted peak. Therefore, one could 
presume that the QND measurements of an arbitrary A^-qubit 
state could be achieved by analyzing the transmission spectra 
of the dispersively-coupled cavity: from the locations of the 
central frequencies of the detected peaks, one can determine 
which basis states {|fc)jv} are superposed; and from the rel- 
ative heights of the corresponding peaks, one can determine 
the superposed probabilities {1/3^1^}. 

Discussions and Conclusions.- Our proposal is based on an 



important assumption, i.e., each detection should be finished 
sufficiently-fast such that the influence of the decay of the de- 
tected quantum state is negligible. This condition is satisfied 
in the current circuit-QED experiment ||9l [lOl uM, wherein 
each data for recording the transmission event of light through 
the cavity can be obtained in about 40 ns and the decay time 
of the detected CPB-qubit is, e.g., Ti ~ 1/is. Thus, the pre- 
dicted multi-peak transmitted spectra could be verified, once 
the CPB-qubits are prepared at the superposition of various 
possible basis states (even the mixed ones). 




FIG. 3: (Color online) Left: (a) Occupation probabilities of various 
basis states versus the decay time for the initial state |11), (b) Spec- 
tral distributions versus the time for the 2-qubit register prepared at 
the basis state |11). Right: (c) Time-averaged transmitted spectra 
over the time r — 0.5 /xs for different initial states: |00) (black), 
jOl) (blue), jlO) (red), and |11) (pink). Here, the decay rate of the 
CPB-qubit is taken as l/Ti,^ ~ 1 MHz. 



Immediately, our proposal could be utilized to explain the 
multi-peak spectra observed in the recent experiment [12]. 
Strictly speaking, the time-averaged (over a relatively-long 
time, i.e., 0.5 /iS, the half of the decay time of the qubit) spec- 
tral distributions shown there are not the desirable readouts 
of the two CPB-qubits, due to the significant decays of the 
detected states during these relatively-long measuring times 
for averages. In Fig. 3(a) we typically show how the prob- 
abilities of various basis states change with the time from 
the decays of the excited state |11) |10|,[lll]. Clearly, dur- 
ing the decay (e.g., at 0.5 /is) the superpositions of the basis 
states are induced. Thus, based on our proposal beyond the 
mean-field approximation, the transmitted spectra would re- 
veal a relevant multi-peak structure. Phenomenally, the decay- 
time dependent steady-state spectrum S*!* {ujl , r) could be ob- 



tained by replacing the unchanged expectable values of the 
qubit-operators in Eq. (8) (i.e^^^(O) and z[^l{0)) with the 
decay-time dependent ones lllC)l[llll . Consequently, the time- 
dependent spectral distributions (due to the decay of the ini- 
tial state 1 11)) was simulated in Fig. 3(b), which agrees basi- 
cally with the corresponding experimental observation (i.e.. 
Fig. 1(E) in Ref. [12]). Thus, except for the non-decayed 
ground state |00) (which coiTesponds certainly to a single 
transmitted peak), the time-evolution spectra due to the de- 
cay of arbitrary excited state, e.g., |10) (or |01), |11)) would 
reveal two (or two, four) peaks. By integrating the decay- 
time dependent steady-state spectrum S'|''(a;L,r), Fig. 3(b) 
shows the relevant time-averaged spectra (over the time inter- 
val T = [0, 0.5] /is) for these time-evolutions. One can see 
that the locations of the averaged peaks agree well with the 
experimental observations 1,12] . While, the relative heights of 
the peaks (marking the basis states induced from the decays 
of the input states) are relatively low. This is an inevitable 
deduction of multiple QND measurements performed sequen- 
tially within the averaged time. 

In summary, an efficient approach to implement the QND 
joint measurements of the N qubits are proposed by detecting 
the transmitted spectra through the dispersively-coupled cav- 
ity. These measurements are the IPMs of the qubits, and thus 
the relevant fidelities could be sufficiently high. Remarkably, 
our proposal is a theory beyond the usual mean-field approx- 
imation and thus the statistical quantum correlations between 
the cavity and qubits are important. In deed, by specifically 
solving the dynamical equation of the cavity-qubit coiTela- 
tions, e.g., for A^ = 1 case, one can prove that the lifetime of 
the qubit-cavity quantum correlation is at the same order of the 
cavity-self. Therefore, the effects of the cavity-qubit coiTela- 
tions, i.e., the transmitted spectra with multiple peaks, could 
be verified by inputting the superposed states of the qubits. 
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